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Abstract 



We study the energy density of two distinct fundamental monopoles in 
SU(3) and Sp(4) theories with an arbitrary mass ratio. Several special limits 
of the general result are checked and verified. Based on the analytic expression 
of energy density the coefficient of the internal part of the moduli space metric 
is computed, which gives it a nice "mechanical" interpretation. We then 
investigate the interaction energy density for both cases. By analyzing the 
contour of the zero interaction energy density we propose a detailed picture 
of what happens when one gets close to the massless limit. The study of the 
interaction energy density also sheds light on the formation of the non-Abelian 
cloud. 
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I. INTRODUCTION 



Since the pioneering work of t'Hooft and Polyakov p] a quarter of a century ago, the 
study of magnetic monopoles in various Yang-Mills-Higgs theories has become a fruitful 
direction of research. Investigations of these solitonic states have uncovered many deep 
and beautiful structures of gauge theories and greatly improved our understanding of those 
theories. 

Magnetic monopoles arise when the Higgs configuration has a non-trivial topology at 
spatial infinity. For a theory with the gauge group G broken into a residue group H, topo- 
logically non-trivial configurations are possible when the second homotopy group of the 
vacuum manifold, namely Il 2 (G/if), is non-trivial. Many works have been performed to 
understand the structure and the metric properties of BPS monopole solutions. It is known 
that when the unbroken gauge group H is Abelian, generalization from the original single 
SU (2) monopole to multi-monopole systems in arbitrary gauge theories is quite straightfor- 
ward (at least conceptually) 0. 

When the unbroken gauge symmetry is non-Abelian, however, the situation becomes 
more complicated. Certain fundamental monopoles (namely the monopoles associated with 
simple roots of the gauge group) become massless and two cases need to be distinguished: the 
total magnetic charge carried by monopole is non-Abelian (as the long range magnetic field 
transforms non-trivially under unbroken symmetry) or purely Abelian. In the former case, as 
was discussed in [||, various topological pathologies appear and prevent us from defining the 
non-Abelian charge globally. On the other hand, when the total magnetic charge is Abelian 
(the latter case), there's no topological obstacle, everything behaves nicely. Therefore the 
majority of the works on monopoles in the presence of non-Abelian unbroken symmetry 
focus this case [|| @ @- The modern picture of such a case is described by the 
so called non-Abelian cloud arising from the interaction between massless monopoles and 
massive monopoles. 

In spite of the progress in understanding the field configurations and the moduli space 
metrics, the detailed behavior of the interaction that accounts for the formation of the non- 
Abelian cloud is still unclear. Under the massless limit, a single monopole will spread out 
and eventually disappear. This trivial behavior can be significantly changed in the presence 
of massive monopoles. In the case where the system carries Abelian charge, the would-be 
massless monopole will lose its identity as an isolated soliton once its core region overlaps 
with massive monopoles, its size will cease to expand, and its internal structure will change in 
a way that reflects the restored non-Abelian symmetry. This picture must be distinguished 
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from the case when the system carries non-Abelian charge (in the massless limit). Having 
a proper detailed description of the two situations will be helpful in understanding the 
formation of the non-Abelian cloud. In this paper we will try to address some of these 
issues. We will compare the behavior of two monopole systems in 577(3) and Sp(4) theories 
since they are the two simplest models containing the interesting contents we are going to 
study. 

The paper is organized as follows: In Sec II, we introduce (as the foundation of our 
calculation) Nahm's formalism for the monopole energy density. In Sec III and IV, the 
energy densities of two (distinct) fundamental monopoles in SU(3) and Sp(4) theories are 
calculated and verified in several special cases. In Sec IV we compute the internal part of the 
moduli space metric from a "mechanical" point of view. In Sec V we study the formation of 
the non-Abelian cloud by analyzing the behavior of the interaction energy density. In Sec 
VI we conclude with some remarks. 



II. NAHM'S FORMALISM FOR THE ENERGY DENSITY 



As was used in many papers, Nahm's formalism has proved to be a powerful tool in cal- 
culating many aspects of monopoles. This method is an analogue of the ADHM construction 
used in instanton physics [|§ and was first proposed by Nahm |H| . Recently Nahm's formal- 
ism has been generalized to deal with calorons (periodic instantons) |l I] [12], and we will 
use some of the results developed in those works. 

Consider the SU (N) Yang-Mills-Higgs system, assuming the asymptotic Higgs field along 
a given direction to be 0^ = diag(/i l5 fi N ) (with //!<•••< /Mr), then the Nahm data for 
the caloron that carries instanton number k are defined over intervals (/ii, fi2),---, (hn-i, Hn), 
(Pn, fJ'i + 2vr//3) (where f3 is the circumference along the S 1 direction in space i? 3 x S 1 ) with 
fii and fii + 2Tr//3 identified. In each interval the Nahm data are triplets of k x k Hermitian 
matrix functions T(t) determined by Nahm's equations and boundary conditions. In this 
paper we will only use the case with k = 1 for which Nahm data are triplets of constants 
representing the positions of the corresponding constituent monopoles. It is known that the 
action density of instantons (in the usual ADHM method) is given by (it differs from [13| 
by a sign since we choose to be Hermitian rather than anti-Hermitian) : 

Ps = trF^ = nniogdet/ (1) 

where / is the inverse operator (whose matrix elements form Green's function) of A^A (A 
is the usual ADHM matrix), and □ is a four dimensional Laplacian. For SU(N) calorons 
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similar results can be established using the Fourier transformation of the original ADHM 
method and one has the following formula for the Green's function f(t,t') = (t\f\f) (11 



^ + |x - T(t)\ 2 +J2 \Ti - T^lSit - tn 



f(t,t')=S(t-t') 



(2) 



where the sum is taken over all the boundary points between adjacent intervals. In order 
to get the Green's function for monopoles instead of calorons, notice that in the constituent 
monopole picture of calorons, an additional type of monopole has been introduced to neu- 
tralize the magnetic charge, so the usual multi-monopole Green's function can be obtained 
by moving the additional monopole to spatial infinity which leads to the following natural 
boundary conditions: 



f(f*i,t') = f(fJ-N,t') = 



(3) 



On the other hand, for purely magnetic configurations, the energy density is given by (a 
factor of 1/2 is omitted for simplicity): 



P = Ps = DD logdet/ 



(4) 



For later convenience, it would be useful to explore Eq.@ a little bit further. Notice 
that for an operator /: 



log det / = tr log / = — tr 
therefore (i = 1,2,3) 

{oo 
EC 1 -/)" 
_n=0 



E 

n=l 



n 



(5) 



dj = tv(r%f) 



(6) 



Further notice that / = (A^A) 1 (Nahm's construction operator A^ is defined to be id t — 
i(x — T) • a in each interval), so 

tr(r 1 9 l /)=tr(AtA^/) 

= tr[^(AtA/)-/9 ?; (AtA)] 

= -E[/<f(*W|x-T,f] (7) 

j 

So finally we have obtained a convenient formula for the energy density: 

dt log det / = - E fe|x - Tj\ 2 J dt'fit', 0] (8) 
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where c^|x — Tj\ 2 has been moved out from the integration since in each interval Tj does 
not depend on i! (for the k = 1 case only). Equation (||) together with Eqs.(|2|) (|]) (f|) 
(in Eq.(^) use three dimensional Laplacian A to replace □) forms a framework to compute 
energy density of monopole systems considered in this paper. 



III. TWO FUNDAMENTAL MONOPOLES IN 5*7(3) THEORY 
The root diagram of 577(3) theory is shown in Fig 1: 
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FIG. 1. Root diagram of SU(3) theory 

We will consider the system formed by one a and one (3 monopole, h in the graph refers 
to the Higgs direction along which f3 monopole is massless. 



A. Energy density 

We choose </>oo (along a given direction) to be diag(— 1 — /a, 2/t, 1 — //) (with — 1/3 < /t < 
1/3), so the masses of the two fundamental monopoles are; 



1 + 3/t; mp = 1 — 3// 



(9) 



Without losing generality, we can place the a-monopole on the origin and /5-monopole on 
(0,0, V) which is equivalent to choosing T(t) = (0,0,0) for t G (—1 — //, 2/t) and T(t) = 
(0, 0, V) for £ e (2//, 1 — //). Applying Eqs.(D) @ to this case we have 



-|^ + |x-T(t)| 2 + X?5(i-2/i) 



f(t,t') = 5(t-t') 



(10) 



- ^ = /(i - = o (ii) 

It's easy to see from these equations that the Green's function has the following form: 
• Case A: -1 — // < t' < 2// 



Asinh[r(t+ 1 + //)] (-1 - /t < t < t') 
/(t, t') = <| B sinh(rt) + C cosh(rt) (t' <t< 2/t) 
D sinh[r'(l -fi-t)] (2fi< t < 1 - fi) 



(12) 



6 



Case B: 2/i < t' < 1 — /i 
f(t,t') = 



A' smh[r(t + 1 + fi)] (-1 - p, < t < 2fi) 

B' sinh(r't) + C cosh(r't) (2// < t < t') (13) 
^ D' sinh[r'(l -//-*)] (f < t < 1 - //) 



where r = y xf + xf, + x|, r' = y xf + a; 2 . + (£3 — ^P) 2 are distances from two monopoles and 
the coefficients A, B, C, A, B', C are all functions of t'. 

In each case Eq. ( PH|) also implies the usual boundary conditions (which we won't bother 
writing down) concerning the continuity of f(t,t') and the jumps of dtf(t,t') at each point 
where the argument of the (^-functions becomes zero. All the coefficients can be computed 
from those boundary conditions. It's helpful to notice that Eq.(||) makes use of f(t,t') only 
in the form of J dt'f(t', t') which is equal to 

/ dt'f(t', t') = [ 2>1 dt'Ait') sinh[r(t / + 1 + //)] 

J J-l-n 

+ f 1 " dt'D'tt') sinh[r'(l - y, - t')] (14) 

J2fi 

and so we only need A and D' . Computing them using boundary conditions and putting 
into Eqs. ([14]) and (H) one obtains the following result: 

rp sinh p sinh q + A\ (p cosh p — sinh p) 



9i(logdet /) = -fi 



rM 



,r'q sinh q sinh p + ^(g cosh q — sinh g) 

" ri r^M (15) 

where p, q, Ax, A 2 , M are defined as: 

p = m Q r; q = mpr' (16) 

Ax = T> sinh q + r' cosh q; A 2 = V sinh p + r cosh p (17) 

M = T> sinh p sinh g + r cosh p sinh q + r' sinh p cosh g (18) 

From Eq. (^) one can also derive the regularized determinant of Green's function to be 

(det/) reg = ^ (19) 

which is defined in the sense that it is finite and gives the same 9j(logdet/) and energy 
density p through: 



7 



^(logdet/) =a i [log(det/) rcg ] (20) 

p = AAlog(det/) reg (21) 

Three typical configurations are shown in Fig 2 (we plot it on the x — z plane since the 
configurations are axially symmetric). 
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FIG. 2. Energy density of two fundamental monopoles on the x — z plane in SU(3) theory, a 
and P monopoles are located at (0,0) and (0,5). The mass ratio m a /mp (from top to bottom) is 
chosen to be 1, 1.35 and 3 
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B. Various limits of the energy density 



In this subsection let's check certain limits of the general form of the energy density, this 
serves as a partial verification of the result obtained in last subsection. 

1. V — case (two monopoles are on top of each other): 

In this case we expect (and will see) that the resulting energy density is the same 
as that of an 5*?7(2)-embedded monopole with mass m = m a + mp = 2. Since two 
monopoles are overlapped, so r = r' and one has: 

M = rsinh(2r) (22) 

So (using Eqs.(|l9D(|0D) 

9i(bgdet/) = -fiH(2r) (23) 

where H(2r) is the m = 2 (m is the mass parameter) case of the single monopole 
function defined as 



Hfmr) = m 



coth(mr) 

mr 



(24) 



From Eq.(p3"D, one can further get: 

A (logdet/) = Wogdet/) = H 2 (2r) - 4 (25) 

and therefore 

p = A A (log det /) = A [H 2 (2r )] (26) 

which is fully compatible (in the suitable convention of normalization) with another 
well known formula fl4] : 



P = A(tr0 2 ) (27) 

since cf) oc H(2r) for single S , f/(2)-embedded monopole with mass m = 2. 
2. Massless limit: 

This is the case when one of the monopoles becomes massless (we will investigate 
this limit in more detail in Sec IV). In our convention this happens when /i = ±1/3. 
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Without losing generality, we choose \i = 1/3 (so the /3-monopole is massless), then 
p = 2r, q = and 

M = r'sinh(2r) (28) 

and therefore 

9 i (bgdet/) = -f i H(2r) (29) 

which (as the T> = case) again leads to the energy density of a single SU (2) monopole 
with mass m — 2. Notice m a = 2 in this case, so such a result means that the j3- 
monopole doesn't affect the energy density in the massless limit (which, as we will see, 
is very different from Sp(A) case). 

3. Removing one monopole: 

Again, without losing generality, let's move the /3-monopole away so that r' — > oo, the 
dominant term in M is 

M ~ r'exp(m / gr') sinh(mQ,r) (30) 

which leads to 

<9j(logdet /) = -fiH(w a r) - m^f • (31) 

Since r[ represents a constant vector at this limit, so Eq.([H]) is exactly what one 
expects for a single monopole with mass m a . 



IV. TWO FUNDAMENTAL MONOPOLES IN 5P(4) THEORY 

Sp(A) (or equivalently SO (5)) theory is the simpest theory to study the non-Abelian 
cloud. The root diagram is shown in Fig. 3: 
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Two cases containing massless monopoles have been studied before. When the asymp- 
totic Higgs field is along the h direction, the Abelian configuration is made of two massive 
/3-monopoles and one massless a-monopole 0, when the asymptotic Higgs field is along the 
h! direction, the Abelian configuration is made of one massive a-monopole and one massless 
/5-monopole |J. In this section we will consider the general energy density for a one a one 
/3 monopole system (with arbitrary mass ratio) which will help us to see how the massless 
limit would be eventually achieved. 



A. Energy density 



We choose 0^ (along a given direction) to be diag(— 1, — fi,fi, 1) (with < \i < 1), so 
the masses of the fundamental monopoles are: 

m a = 2 — 2/i; = 2/x (32) 

The monopole locations are chosen to be the same as in the SU (3) case. As we know, 
in Nahm's method Sp(4) theory is embedded into SU(4) theory whose Nahm data satisfy 
symmetric constraints 0, in our convention the SU{4) Nahm data are given by T(£) = 
(0, 0, 0) for t E (-1, -fi) U (/i, 1) and T(t) = (0, 0, V) for t e {-fi, //). Applying Eqs.(|) © 
to this case we have: 



~ + |x - T(t)| 2 + V5(t + fi) + V5(t - fi) 



f(t,t') = 5(t-t') 



/(-i,0 = /(i,0 = o 

The Green's function satisfying these equations has the following form: 
• Case A: -1 < t' < —fj, 

' Asinh[r(t + 1)] (-1 < t < f) 

B sinh(rt) + C cosh(rt) (t' < t < -/x) 
Dsmh(r't) + E cosh(r't) (-/i < t < /i) 
t Fsinh[r(l -*)] (ji<t<l) 



Case B: —fj, < t' < fi 



f&t) 



A'sinh[r(t + 1)] (-1 < t < -//) 

B' sinh(r't) + C cosh(r't) (-// < t < t') 

D' sinh(r't) + E' cosh(r't) (f < t < /j) 

F'smh[r(l-t)) (ji < t < 1) 



(33) 



(34) 



(35) 



(36) 
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• Case C: // < if < 1 



Z(M') 



A" sinh[r(i + 1)] (-1 < t < -//) 

5" sinh(r't) + C" cosh(r't) (-/i<t< fi) 
D" sinh(rt) + cosh(rt) (n < t < t') 
^ F" sinh[r(l - t)] {if < t < 1) 



(37) 



where r and r' have the same meaning as before, and A, F" are all functions of t'. 

Similarly as in the SU(3) case, we only need some of the coefficients (A, B', C, D") 
which (as well as all other coefficients) can be determined by boundary conditions coming 
from the 5-functions in Eqs.(|3"3|). Compute them and putting them into Eqs.(|8|) one obtains: 

LB 2 



di (log det y) = — T{ ~ — ~~ — ~ — f*' { — ~ — ~~ 

* cosh rP + sinh rQ 1 \ L 

where the following functions are introduced: 

u = (1 — /i)r; v = fir; w = fir' 



M 



(38) 



(39) 



— coshr + cosh(w — v) , . 

Ai = i + 2 1 - //) sinhr 

r 

— sinhr — sinhf-u — v) , . 

A 2 = + 2(1 - ji) coshr 



(40) 
(41) 



sinh(2 / ur / ) sinh(2/ir') 
B 1 =fi-\ — ; B 2 = fi- 



2r> 



2r' 



(42) 



L = T> sinh u cosh w + r cosh u cosh w + r sinh u sinh w 
M = T> sinh u sinh w + r cosh w sinh + r' sinh u cosh w 
Ni — T> sinh u sinh w — r cosh u sinh w + r' sinh i> cosh w 
N 2 = T> cosh f sinh w — r sinh t> sinh w + r' cosh t> cosh w 
N s = T> sinh u sinh(2w) + r cosh u sinh(2w) + r' sinh u cosh(2u;) 



(43) 
(44) 
(45) 
(46) 
(47) 



P = r' sinh vM - iVi JV 3 ; Q = -r' cosh vM + N 2 N S 



The reguralized determinant of The Green's function from Eq. 

(det /) reg 



is 



2 / 



LM 



Three typical configurations are shown in Fig 4. 



(48) 



(49) 
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FIG. 4. Energy density of two fundamental monopoles on the x — z plane in Sp(4) theory a 
and (5 monopoles are located at (0, 0) and (0, 5). Mass ratio m a /mp (from top to bottom) is chosen 
to be 1, 2 and 4 

One can see that the energy density of the two fundamental monopoles in 5p(4) theory is 
not symmetric, this is because here a and (5 monopoles have different energy distributions. 
As we know, the mass of SU (2)-embedded 7-monopoles in arbitrary gauge theory is given 
by h ■ 7* while the scale of such monopoles is determined by l/(h • 7) (7 is the root to which 
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the S77(2)-embedded monopoles are associated, h is the asymptotic Higgs direction). As a 
consequence, in Sp{4) theory when a and (3 monopoles have the same mass, the scale of the 
/3-monopole is only half of the scale of the a-monopole therefore the maximal energy density 
is eight times that of the a-monopole, this is why in the plotting with m a : mp = 1 one can 
hardly see the a-monopole. 



B. Various limits of the energy density 

The expression of the energy density of the two fundamental monopoles in Sp{4) theory 
is much more complicated than in the SU(3) case, to see its correctness, let's check several 
special cases. 

1. T> — case (two monopoles are on top of each other): 

In this case we again expect an S"?7(2)-embedded monopole with mass m = m a + mp = 
2. Since r = r', one gets 

L = rcoshr; M = rsinhr (50) 



which generates the familar result: 

9i0ogdet/) = -r i H(2r) (51) 



which is the same as Eq. (|2"3"D , so in this limit the result is exactly what we expected. 
2. Massless limit: 

We are interested in the limit when the total magnetic charge is Abelian which is 
realized by [i — > (when the /3-monopole becomes massless). In this limit one has 

L = T> sinh r + r cosh r; M = r' sinh r (52) 

which leads to 

, „ N 2r 2 cosh(2r) + 2Vr sinh(2r) - r sinh(2r) - 2V cosh(2r) + 2V . . 

aj(logdet /) = — Ti — — r — — — — r — — (53) 

H sinh(2r) + T>r cosh(2r) — T>r 

On the other hand, the Higgs configuration in this limit is given in [0 ||, in that 
work, the Higgs configuration of N — 1 fundamental monopoles is calculated in the 
symmetry breaking pattern SU(N) — > U(l) X SU(N - 2) X U(l). The relevant Sp(A) 
Higgs configuration can be obtained by taking the N = A case and putting two massive 
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monopoles on top of each other. In a proper normalization the result of 
simplification) can be written as (where a r = a ■ f) 



1 

71 



H(2r)a r 



2T> tanh 2 r 



2£>tanh 2 r 



sinh(2r) (r+X> tanh r) 
2Z>r-£>sinh(2r) 



y V sinh(2r)(r+Dtanhr) 2r cosh 2 r(r+V tanh r) °"' r J 



(after 



(54) 



which leads to 



tr0 2 = H 2 (2r) + - 



AV tanh 2 r 



sinh(2r)(r + Dtanhr) 



+ 



2Vr - £>sinh(2r) 
2r cosh 2 r (r + V tanh r) 



(55) 



As we did before, in order to show that Eqs.(§J) (|53D and ( p7|) ( |55|) give the same result 
it is sufficient to show that the difference between <9i<9i(logdet /) and tr0 2 is constant. 
This is messy, but one can verify that 



(9,(9, (log det /) =tr0 2 -4 



(56) 



so the massless limit works out correctly. Unlike the SU(3) case, the massless limit of 
Sp(4) theory is rather non-trivial (the non-Abelian cloud coming into playing). 

One can easily show that when the a-monopole becomes massless (in that case 
the total magnetic charge is non-Abelian), the situation is similar to the SU(3) case, 
namely the energy density is equal to the energy density of a single /3-monopole. 

3. Removing the a-monopole: 

Now let's check what happens when one monopole is removed. Since the two monopoles 
are not symmetric in the Sp(A) case so we check them separately. When T> — » oo 
but keeping r' finite (therefore r — > oo), the a-monopole is removed. The leading 
contributions of L, M are: 

L ~ r exp[(l — /j,)r] cosh(//r'); M ~ r exp[(l — /i)r] sinh(/xr') (57) 

Putting these into Eqs.(^) and (|20| ) one obtains: 

di (log det /) = -r'^fir') - m a f< (58) 

This (notice that fj is a constant vector at this limit) represents an S"[/(2)-embedded 
monopole with mass m = which is what we are expecting. 
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4. Removing the /3-monopole: 

When T> — > oo, r' — > oo, the /3 monopole is removed. The leading contributions of 
L, M are 

L ~ M ~ r' exp(//r') sinh[(l - //)r] (59) 

which lead to (again, the second term is a constant vector at this limit) 

^(logdet/) = -2rH[(l - /u)r] - (60) 

The energy density coming from this expression is the same as two directly superposed 
SU(2) monopoles with mass m = 1 — fi = m a /2. This is consistent with the 57/(4) 
picture introduced in the discussion of the massless limit, namely the a monopole 
can be considered as the two overlapped SU (4) monopoles whose energy densities are 
simply added since they are non-interacting. This also gives a direct demonstration of 
our discussion at the end of subsection A. 



C. From the moment of inertia to the moduli space metric 



Since we have the analytic form of energy density, we are now able to compute the 
internal part of the moduli space metric using a nice "mechanical" interpretation. The idea 
of using a mechanical interpretation can be traced to Manton's original work |15[ (where 
the concept of the moduli space metric itself was introduced by comparing the action of 
the monopole system and mechanical system) and was used in certain arguments in recent 



works |6j [16]. In this subsection we will consider the massless limit of our Sp(A) system. 



The metric in this case is known to be || (changed into our convention): 

167T 2 47T 

ds 2 = mdx 2 + —^- d X 2 + -^ dT>2 + ^V(al + a 2 , + a%) (61) 

where m is the total mass (which is just the mass of the a-monopole in this case) of the 
system, <7i,<t 2 and <r 3 are 1-forms defined as 

<j\ — — sin ipd9 + cos ip sin 8d(f) (62) 
cr 2 = cos ipdO + sin ip sin 6d(p (63) 
03 = dijj + cos 9dcj) (64) 

with the Euler angles 9,4> and ip having periodicities tt, 2ir and 4n, respectively. It's in- 
teresting to notice that the last term in Eq.flBTD has a "mechanical interpretation" as the 
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rotational energy associated with the massless cloud with the coefficient AnT> playing the 
role of the moment of inertia of the cloud. To see this let's calculate the moment of inertia 
of the non-Abelian cloud. As we know, for a spherically symmetric system the moment of 



inertia tensor has the form 7jj = I5ij with 

i = l 

3 



J dVr 2 p (65) 



Since the rotation of the cloud is actually the gauge rotation in internal space, we should 
remove the gauge invariant part p(D = 0) (this is an S"?7(2)-embedded 7-monopole which is 
gauge invariant), so the effective energy density relevant to the internal rotation is p(T>) — 
p(0). Using the result obtained in the last subsection one can verify that 

2 



I = - J dVr 2 [p(V) - p(0)] 
drr A [p(V) - p(0)] 



_ 87r 

~ T 

= WttV (66) 

which leads to a term lQ7cT>(dujf + dio\ + du%) in the moduli space metric. Here 1-forms 
dui,du 2 and du 3 are defined in the group space of SU{2) (namely S 3 ). To compare with 
the last term of Eq.(|61|) one notices that 

o\ + a\ + al = d6 2 + dcj) 2 + dip 2 + 2 cos 6d0dip (67) 

therefore the volume of (6,<f),if)) space is (g = det(py-) is the determinant of metric matrix 
coming from Eq.(|67D) 

.2 



V = / ^\g\d6d<pdi) = 16tt 2 (68) 

Since the volume of group space S 3 is 2ir 2 the two sets of 1-forms are related by crj = 2duji 
(i = 1,2,3). Taking this into account we see that lQirV(duj 2 + dio\ + du 2 ) computed from 
the moment of inertia is in accordance with the last term in Eq.flBTD obtained using other 
methods. 



V. INTERACTION ENERGY DENSITY AND THE FORMATION OF THE 

NON-ABELIAN CLOUD 

In previous sections we have calculated the energy density of two monopole systems in 
SU(3) and Sp(4) theory. We've already known that when one approaches the massless limit, 
the situations are very different depending on the total magnetic charge. If the total magnetic 
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charge is non-Abelian, the resulting energy density is simply the same as the energy density 
of the massive monopole. When the total magnetic charge is purely Abelian, however, the 
energy density distribution is deeply affected by the existence of massless monopoles. In the 
latter case, there is a non-Abelian cloud surrounding the massive monopole, neutralizing the 
non-Abelian components of the magnetic charge. In this section we want to have a close 
look at the evolution of the energy density when one approaches the massless limit. 

There's no unique choice of quantity to describe the formation of the non-Abelian cloud, 
nor is there any unambiguous definition of the non-Abelian cloud itself. But physically 
there is no doubt that it is the interaction between massive and massless monopoles that 
determines the behavior of the system, including the formation of the cloud. Our strategy 
is to study the interaction energy density defined as 



where p a and pp are the energy density of isolated a and (3 monopoles. p mt describes the 
change of energy distribution caused by the interaction between two monopoles. In particular 
we will look at the contour of the zero interaction energy density which gives information 
on where interaction gathers energy and from where it extracts energy. 

We have used Maple to generate numerical data and plotted several typical contours 
(shown in Figs 5, 6) for SU (3) and Sp(A) theories (the region enclosed by the contour has 
a positive interaction energy density). 



Pint — Ptotal — Pa — P(3 



(69) 
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FIG. 5. Contour diagrams of the zero interaction energy density on the x — z plane for SU(3) 
theory, a and (3 monopoles are located at (0,0) and (0, 10), mass ratios m a : mp are choosen to 
be (from top to bottom) 4, 7, 19, 199 
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z 

FIG. 6. Contour diagrams of the zero interaction energy density on the x — z plane for Sp(4) 
theory, a and (3 monopoles are located at (0,0) and (0, 10), mass ratios m a : rap are choosen to 
be (from top to bottom) 4, 19, 66, 99 

From the contour diagrams one can see the major difference between the two theories 
when approaching the massless limit. In both cases we start with a simply connected contour 
for the small mass ratio (when the distance between two monopoles is large, the starting 
contour could be different). The contour deforms and grows when the mass ratio increases, 
in both cases it breaks into two disjoint pieces when mass ratio is sufficiently large. The 
reason it breaks can be understood by directly analyzing the massless limit of p int (in the 
SU(3) case p nit itself vanishes but one can use Pint/^/3 which remains finite), Figs 7, 8 show 
those limits: 




-f 1 1 1 i 1 1 1 1 t 



FIG. 7. Pmt/m/3 in SU(3) theory at the massless limit (this curve is independent of V) 
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FIG. 8. pi n t in Sp(4) theory at the massless limit for V = 10 (this curve has a weak dependence 
on V) 

In both cases the limiting interaction energy density becomes negative outside the core of the 
a-monopole, this means that the contour can't keep growing and remains simply connected 
all the way as one increases the mass ratio. After the breaking point, the part of the contour 
that surrounds the ct-monopole is stabilized, but the other part undergoes a very different 
evolution in the two cases: In the SU(3) case that part of the contour keeps growing, but 
gradually moves away from the a-monopole (because of the scale of the graphs, it might not 
be able to see that easily, but the shortest distance between the two parts of the contour 
is increasing). In the Sp(4) case, however, the situation is the opposite, the other part 
of the contour eventually shrinks and finally disappears (this happens before going to the 
limit). It should be mentioned that in Sp(A) case if the /5-monopole is inside the core of the 
a-monopole, the contour could shrink and be stabilized without breaking into two pieces at 
first. 

The difference in the contour diagram of the two cases explains their different massless 
limits and accounts for the formation of the non-Abelian cloud. Although in both cases, the 
interaction alters the energy distribution by accumulating energy in certain regions, in the 
case of SU (3) that region is ever expanding and the effect of the interaction (therefore the 
massless monopole itself) is smeared out over the infinitely large area, so the final energy 
density is completely dominated by the remaining massive monopole. On the other hand, 
in the Sp(4) case, the interaction extracts energy and deposits it into a small region (in 
some sense one can say that the interaction is more "localized" in this result 
it affects the energy density distribution significantly, and because of the interaction the 
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massless monopole doesn't grow into infinite size as it would if isolated. The non-Abelian 
cloud is just the effect of such an interaction. 

VI. CONCLUSIONS 

So far we have analyzed the energy density of two monopole systems in SU(3) and Sp(4) 
theories and obtained some idea on how the massless cloud forms. Based on these results 
one can make some qualitative conjectures on what might happen in general cases where 
the interaction energy can be defined as 

Pint Ptotal Pmassive Pmassless 

(70) 

The basic property of pi nt we learned from the previous observation is that when the total 
magnetic charge is purely Abelian, the interaction is more "localized" in contrast to the 
opposite case. Such an interaction extracts energy from distant regions, accumulates it in 
the vicinity of massive monopoles and gradually builds up the structure of the non-Abelian 
cloud. This is fairly similar to the Sp(4) case. 

When the total magnetic charge is non-Abelian, however, qualitatively different situa- 
tions could arise in general. To see this, let's look at the case with the symmetry breaking 
pattern SU(N) -> £7(1) x SU(N-2) x £7(1) (N > 4). Let au...,a N -i denote simple roots. 
When the system contains massive ai, aN-i and would be massless aij (i = 2, ...,N — 3) 
monopoles (notice that the a7v_ 2 -monopole is absent and so the total magnetic charge of the 
system is non-Abelian), only massive monopoles survive at the massless limit. This can be 
seen by noticing that the system under consideration is equivalent to the system studied in 
H (which contains the ajv-2-monopole as well and so the total magnetic charge is Abelian) 
with the a7v-2- mc, nopole removed. From || we know that the only cloud parameter of the 
system is given by 

N-l 

P=^|x i -x i _ 1 | (71) 

i=2 

So removing any massless monopole is equivalent to removing the whole cloud (since it 
makes the cloud size infinity) therefore only massive monopoles survive. This situation is 
similar to the SU(3) case we have studied. But there are other systems which don't show 
such a direct analogue. As an example we can go back to our Sp(4) theory, and consider a 
system with N massive a-monopoles and N — 1 would be massless /5-monopoles (so the total 
magnetic charge is non-Abelian). At massless limit, the massless monopoles in this system 
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will form a non-Abelian cloud rather than disappearing. This is because such system can be 
obtained by removing one /3-monopole from a system containing iV (N > 1) a — (5 pairs. At 
the massless limit the latter system contains a non-Abelian cloud with many independent 
size parameters. Removing one /5-monopole will not make all these parameters infinity and 
therefore will not destroy the whole cloud. Another way to understand this is to notice that 
when the extra a-monopole is removed, we are left with a system made of N — 1 a — (3 pairs 
(it can be called an Abelian sub-system) which certainly contains the non-Abelian cloud. 
Since removing the a-monopole won't create cloud, the cloud must exist in the original 
system. This argument can be generalized. 

These examples reveal the complexity of the general cases. It seems that at the massless 
limit a system with a non-Abelian total magnetic charge can still contain massless monopoles 
(in the form of a non-Abelian cloud) in a "maximal Abelian sub-system" . We think further 
considerations on such situations will be interesting. 



Note added: While writing this paper, we noticed the appearance of [|17| from which the 
energy density of SU (3) and Sp(4) monopoles can also be obtained. 
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